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Abstract 

The one-dimensional problem of N particles with contact interaction in the presence of a 
tunable transmitting and reflecting impurity is investigated along the lines of the coordinate 
Bethe ansatz. As a result, the system is shown to be exactly solvable by determining the 
eigenfunctions and the energy spectrum. The latter is given by the solutions of the Bethe 
ansatz equations which we establish for different boundary conditions in the presence of 
the impurity. These impurity Bethe equations contain as special cases well-known Bethe 
equations for systems on the half-line. We briefly study them on their own through the 
toy-examples of one and two particles. It turns out that the impurity can be tuned to lift 
degeneracies in the energies and can create bound states when it is sufficiently attractive. 
The example of an impurity sitting at the center of a box and breaking parity invariance 
shows that such an impurity can be used to confine asymmetrically a stationary state. This 
could have interesting applications in condensed matter physics. 
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Introduction 



Forty years ago, E. Lieb and W. Liniger published their seminal paper presenting exact results 
for the one-dimensional repulsive Bose gas [I], extending the previous investigation for hard-core 
bosons [2]. This was completed in [3] for the attractive interaction. It is remarkable that this 
purely theoretical work finds a huge amount of applications nowadays with the advent of optical 
lattices. The latter allow to produce quasi one-dimensional environment where the quantum 
behaviour of ultracold atoms can be probed experimentally [4j. The main ingredient used in 
[1] is the celebrated Bethe ansatz [5j for the wavefunction. In essence, this ansatz assumes an 
expansion of the wavefunction on plane waves and the coefficients are determined so as to take 
the interactions into account. Then, the energy spectrum is given by the solution of the Bethe 
ansatz equations. Soon after, M. Gaudin [6] and then C. N. Yang [7] generalized the results for 
particles with different type of statistics by considering a wavefunction in different irreducible 
representation of the permutation group. In particular, the investigation of C. N. Yang relied 
on the now famous Yang-Baxter equation [8]. Finally, M. Gaudin studied also the analog of the 
system of p] when the bosons are enclosed in a box [9]. In particular, he introduced a slightly more 
general Hamiltonian than the contact interaction Hamiltonian of [T] depending on two different 
coupling constants. The latter was recovered recently in [10], for particles with arbitrary spin, as 
a limit of a long range interacting Hamiltonian of Sutherland type [UJ for which integrability was 
proved. It was also shown that the symmetry of this system is the reflection algebra symmetry 
[T21 [13] . This motivates the interpretation of the Hamiltonian considered by Gaudin as describing 
particles on the half-line, or equivalently in the presence of a purely reflecting impurity. 

Let us stress that the many-body Hamiltonian of [1] is the restriction to the iV-particle Fock 
space of the well-known nonlinear Schrodinger (NLS) Hamiltonian (see e.g. [14] for a review). The 
NLS model is one of most studied examples of integrable field theory for which a huge amount 
of exact results is known. In the same way, the Hamiltonian of [9] is the counterpart of the 
NLS model on the half-line whose symmetry is given by the reflection algebra ^S\, showing the 
consistency of the approach of [10J. In [15], the concept of boundary algebra [TJ] was crucial to 
establish all the properties of NLS on the half-line as an integrable system. 

The question of a reflecting and transmitting impurity in integrable systems appeared naturally 
as a generalization of a purely reflecting boundary. The first approach in this context was done in 
[IT] where a set of reflection-transmission equations was derived. Later in [18], it was proved that 
non-trivial solutions for the two-body scattering matrix do not exist if we require a non-vanishing 
reflection and transmission coefficients for the impurity. More recently, a new framework was 
introduced to handle reflecting and transmitting impurities in integrable systems [19] and was 
shown to be more general [20J. It was successfully applied to the NLS model with impurity 
[2T1 [22] providing the first non-trivial integrable system with impurity. 

Consequently, it seemed natural to us to consider the many-body analog of NLS with impurity 
and to investigate it along the lines of [Tj and [9]. Just like the system without impurity, it may 
be of particular interest for current experiments in condensed matter physics. 

After presenting the problem in Section 1 together with some notations to describe it, we 
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show in Section 2 that it is exactly solvable thanks to an appropriate Bethe ansatz for the N- 
particle wavefunction. In Section 3, the full use of the Bethe ansatz combined with the physical 
requirement of a finite size system allows to establish the Bethe ansatz equations in the presence of 
an impurity. This, in turn, is well-known to determine the energy spectrum. Section 4 is devoted 
to specific examples. First, we show that our setup reproduces the results of [5] as a special case. 
Then, we use the one and two-particle toy examples to illustrate the effects of the impurity 

on the energy levels and on the parity symmetry. Finally, in Section 5, we present our conclusions 
for this work and give an outlook of future investigations. 



1 The nature of the problem 

1.1 Combining two systems 

In this paper, we study a one- dimensional system of N particles interacting through a repulsive 8 
potential in the presence of an impurity sitting at the origin and described by a point-like external 
potential. This problem is the combination of the interacting system studied in [TJ [7] and the free 
problem in the presence of a point-like potential, see e.g. [23] and [23]. 

Each of these problems has a well-defined translation in terms of a partial differential equation 
problem together with boundary conditions for the wavefunction. For example, let us denote by 
<f(xx, . . . , Xjv) the iV-particle wavefunction for a gas with a repulsive 5 interaction of coupling 
constant g > 0. Then, following [lj, cp is solution of the free problem for the energy E 

N 

f(xi,...,x N ) = E <p(x u ...,x N ) , (1.1) 

i=l 

with the additional requirement of continuity and jump in the derivative at each hyperplane 

Xj = Xk, j ^ k 

ip(x l7 ...,x N )\ x . =x + = ip(x u ...,x N )\ Xj=x - (1.2) 

(d x -d Xk ) cp(x u ...,x N )\ + = [(d x -d Xh )+2g]<p(x 1 ,...,x N )\ - (1.3) 



Now, in [23], the second problem is presented for the one-particle wavefunction tp(x), x^O using 
a unitary matrix U G (7(2) characterizing the impuritjQ: 

lim ((U - I)$(x) + i(U + lW(x)) = , (1.4) 

where 

*w=U-t)) • *'w=(_$i)) .*><>. (") 

<p'(x) = d/dx (p(x) and I is the 2x2 unit matrix. The matrix U can be parametrized as follows 
H ( V " 



U = e*[ _^ ) , f G [0,7r), fi,u G C such that \fi\ 2 + |z/| 2 = 1 . (1.6) 



Hn [23j . there is a length scale Lq which is shown to be an irrelevant parameter. We set it to 1 in this paper. 
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The symbol * stands for complex conjugation. Mathematically, this problem corresponds to all 
the possible self-adjoint extensions of the free Hamiltonian when the point x = is removed from 
the line. 

As announced in the introduction, the goal of this paper is to present and solve the quantum 
iV-body problem combining these two models. Physically speaking, we address the problem of 
a one- dimensional gas of interacting particles in the presence of a tunable impurity in the sense 
that the parameters in (11.61) are free and can therefore be used to model different impurities with 
different coupling constants. 

1.2 Notations and definitions 

From the mathematical point of view, the lesson we learn from [U [7] is the crucial role played by 
the permutation group &n of N\ elements. It consists of N generators: the identity Id and N — 1 
elements Ti, . . . , Tat_i satisfying 

TjT 5 = Id , TjTe = T(Tj for \j-£\>l , (1.7) 

In particular, the last relation gives rise to the famous Yang-Baxter equation [7J, [8]. For conve- 
nience, we denote a general transposition of 6^ by T^-, i < j, given by 

Tij = Tj-i ■ ■ ■ Tj + iTjTj + i . . . Tj-i (1-9) 

Then, in [9], the role of the so-called reflection group was emphasized and in [16], the Weyl 
group 2Ujv associated to the Lie algebra B N replaced the permutation group in the construction 
of a Fock space for systems on the half-line. Let us note that the same group proved to be 
fundamental in the constructions of [25] corresponding to an interacting gas on the half-line where 
the usual 5 interaction was replaced by another contact interaction, the so-called 5' interaction. 
2Uat contains 2^^! elements generated by Id, Ti, . . . , Tjv-i and R\ satisfying ( 11.71) . ( 11.81) and 

R 1 R 1 = Id , (1.10) 
R 1 T 1 R 1 T 1 =T 1 R 1 TiR 1 , (1.11) 
R 1 T j = T j R 1 for j > 1 . (1.12) 

Let us define also Rj, j = 2, . . . , N as 

R j = T j _ l ...T 1 R 1 T 1 ...T j _ 1 (1.13) 

Remarkably enough, the same group appears in the construction of Fock space representations 
for systems with an impurity in the context of RT algebras [19] . One may wonder how the same 
structure can account for systems on the half-line (i.e. with purely reflecting impurity) and also for 
systems on the whole line with a reflecting and transmitting impurity. The essential point is the 
choice of representation. Typically, for a system on the half-line involving particles with n internal 



3 



degrees of freedom, n-dimensional representations of Wn are used. It was realized in [211 122] that 
the same problem on the whole line with a reflecting and transmitting impurity requires 2n- 
dimensional representations of 2Uat. The interpretation of this fact is that the impurity naturally 
defines two half-lines which are physically inequivalent. Thus, in addition to the n degrees of 
freedom of the internal symmetry, each particle carries two degrees of freedom + or — labelling 
the side of the impurity. These correspond to the two components of the wavefunction $ in (12. 2p . 



2 Exact solvability of the model 

For pedagogical reasons, we present first the one and two-particle cases in detail before turning to 
the study of the iV-particle problem in its full generality. We refer the experienced reader directly 
to section 12.31 

2.1 One particle 

For i6l \ {0}, the one-particle wavefunction is taken as follows 

<PW = < _. , (2.1) 
I tp \x) x < (J 

Note that no parity relation is assumed for tp, which is a crucial point in our approach. We define 
for x > 

•<*> - ( ) <«> 

and following the previous paragraph, the boundary conditions at x — 0, which we will call in this 
paper the impurity conditions, read 

(U-l)$(x)=-i(U + l)&(x) for x^0 + , (2.3) 

the matrix U being given in (11.61) . Let us expand $ on plane waves as follows 

$(x) = exp(ikx)Aid + exp(— ikx)Au (2.4) 

where Ap = ( j^L J, P = Id,R. These coefficients are constrained by condition (12.31) . This is 

essentially the celebrated Bethe ansatz for one particle and it is solution of equation (11.11) with 
E = k 2 . Plugging back into ( 12. 31) . one gets, 

A R = Z(-k)A Id and A Id = Z(k)A R (2.5) 

where 

Z(k) = —[U — I — k(U + I)}' 1 [U - I + k(U + I)} (2.6) 



4 



The consistency of the ansatz is ensured by Z(k)Z(—k) = I which is readily seen to hold. The 
property Z'(k) = Z(—k), where > stands for Hermitian conjugation, then leads to the physical 
unitarity Z\k)Z{k) = I. 

For completeness, let us make the connection with the other usual setting of the problem. For 
v 7^ 0, (12.31) is equivalent to 

<p(x) \ ( a b \ / <f(— x) 



y(x) j- a {c d JU'(-4 l> fOT ^ 0+ > < 2 - 7 > 

where 

{a, d e K, a e C : ad - be = 1, a*a = 1} . (2.8) 

This is the SU{2) parametrization. Writing n = ji R + ifij, v = v R + ivj with fi R , /i/, is R , ujE R, 
the relation between the two parametrizations is 

zV _ sin£-/jj ^_ cos£ + /i fi _ cos£ - //g rf _ sin g + /x f 

|z/| |z/| 





, a = 











From this one finds 



T-(-jfe) fl-(-ife) 



where 



6fc 2 + z(a-d)fc + c + 2za/c f2 nl 

1 J + z'(a + d)Jfe - c ' 1 J M; 2 + z(a + d)£;-c' K ] 

_ bk 2 + t(a-d)k + c _ -2ia*k 

1 ' bk^-i{a + d)k-c 1 1 ' bk^-i{a + d)k-c 1 { ' 

are usually referred to as reflection and transmission coefficients of the impurity. Of great impor- 
tance is the well-known associated basis of orthonormal eigenfunctions for scattering states 

ij,+(x) = 6{-x)T-{k)e lkx + 6{x) [e ikx + R + {-k)e~ ikx ] , k < , (2.13) 
tJ,-(x) = 6{x)T + {k)e ikx + 0(-x) [e ikx + R~ \-k)e' ikx ] , k > , (2.14) 

which appears as a particular choice of the above setting, justifying the Bethe ansatz approach. 
These eigenfunctions play a crucial role in the quantum field theoretic version of this problem i.e. 
the nonlinear Schrodinger equation with impurity |21j . 

For v = 0, (12.31) gives rise to the so-called separated boundary conditions of the form 

^+'(0+)=g + ^ + (0+) , <p-'(-0 + ) = q- <p-(-0 + ) (2.15) 

with q + , q~£ R U {oo} given by q* 1 = =Ftan(^0, C being the argument of fi. 
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ip(x 1 ,x 2 ) 



(2.16) 



2.2 Two particles 

In the same spirit as before, for x±,x 2 6 R \ {0} and x\ ^ x 2 , the two-particle wavefunction is 
taken to be 

' ip ++ (xi, x 2 ) X\ > 0,x 2 > 
ip + ~(x 1 ,x 2 ) x 1 >0,x 2 <0 
ip^ + {x\ 1 x 2 ) x\ < 0,x 2 > 

{(p (xi,X 2 ) X!<0,X 2 <0 

Then, we define for x±, x 2 > and x\ ^ x 2 

( (p ++ (x 1 ,x 2 ) \ 
(p + ~(x 1 , -x 2 ) 

(p' + (-X 1 ,X 2 ) 

\ y—{-x 1 ,-x 2 ) I 

Now, we implement the fact that each particle can interact with the impurity by imposing two 
impurity conditions 

[{U -l)®l]^{x u x 2 ) = -i[(U + l)®l]d xi ${x 1 ,x 2 ) for x ± -> + (2.18) 
\1®(U -T)]$(x 1 ,x 2 ) = -i\l®(U + T)]d xa $(x 1 ,X2) for x 2 -> 0+ . (2.19) 

The interaction in the bulk between the two particles through a 5 potential is implemented as 
follows 



^(x 1 ,x 2 ) 



(2.17) 



$ (^2)U = 4 = T x $( Xl ,x 2 )\ xi=x - 
(d Xl -d X2 )<$>( Xl ,x 2 )\ Xi=4 = f 1 [(d xl -d X2 ) + 2g]^(x 1 ,x 2 )\ xl=x - 

where Ti is the representation on C 2 (g) C 2 of Ti E & 2 given by 

/ 1 \ 

T\ = 



(2.20) 
(2.21) 



10 
10 
\ 1 ) 



(2.22) 



Similarly, Id is the 4x4 unit matrix representing Id. 

The crucial and new point now is to formulate an ansatz for <&(x\,x 2 ) and show that it solves 
the problem. For < xq\ < xq 2 with Q <E & 2 = {Id, Ti}, we take 

$q(x 1 ,x 2 )= ^2 exp[i(k P1 x Q1 + k P2 x Q2 )] QA P {Q) (2.23) 



Pe w 2 



( A+ + {Q) \ 
A + p~{Q) 
A~ P + (Q) 

\ A- P ~(Q) I 

The impurity conditions imply 



where A P (Q) 



are the coefficients to determine. The energy is simply E = kj + k^. 



A PRl (Id) = [Z(—kpi) ® I] Ap(Id) 
T\ Apr^Ti) = [I®Z(-k P1 )\ fx A P (T X ) 



(2.24) 
(2.25) 
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which reduce to 



A PRl (Q) = [Z(-k P1 ) <g> I] A P (Q) with Q e 6 2 



using 



Ti [I®Z(fc)] Ti = Z(&)®I. 

The matrix Z is the one given in (12.61) . 

The bulk conditions fl2T20l) and Q72TT > give 

1 



fepi - ^P2 + ig 
Introducing the eight-component vector 



Ap(Id) 

A P (T!) 



We can rewrite (I2.26P and (12.281) in a compact form 

Apr, = [I® Z(-fc P1 ) 
*4 PTl = F(/c P1 - 

where 



51] .Ap 



k+ig 



k+ig 



-if I 



k+ig 



k+ig 



k+ig 



k+ig 



\ 



k+ig 



k+ig 



\ 



k+ig 



k+ig 



k+ig 



k+ig 



-w 



k+ig 



-if) 



k+ig 



-ig 



k+ig 



J 



(2.26) 



(2.27) 



{{kpi - k P2 )Ap(QT 1 ) - igA P {Q)) with Qe 6 2 (2.28) 



(2.29) 



(2.30) 
(2.31) 



^2.32) 



Since the relations R\ = Id, 7\ 2 = Id and R%TiRiTi = TiR{TxRi hold in 2U 2 , equations (12.301) 
and (I2.3ip require that Y(k) and Z{k) satisfy the consistency relations 

Z{k)Z{—k) = I , Y(fa - k 2 )Y(k 2 - fa) = I® I® I (2.33) 

and a generalization of the celebrated reflection equation p2J HB], 

Y(u - v)[I <g> Z(u) ® I]y(w + t))[I® Z(v) ® I] = [I® Z(u) ® l]Y(u + v)[I® Z{u) ® I]F(?i - u) 

The explicit form of F and Z ensures the validity of these equations. 

It is a generalization in the sense that even in the scalar case (particles with no internal 
degrees of freedom), our setup produces a two-dimensional representation of W 2 . This is the first 
illustration of the general statement at the end of Section [U 

We conclude that the two-particle model is exactly solvable in the sense that the eigenfunction 
can be consistently given starting from a given Ap, say Am- 
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2.3 N particles 

Following the previous arguments, we present the general solution of our problem for N particles 
and prove its exact solvability. 

For xx, . . . ,xn G R \ {0} and xi, . . . , x^ 2 by 2 different, the natural generalization of (12.161) 
for the wavefunction is 

ip(xi, . . . , xjy) = ip ei '" €N (xi, . . . , xn) in the region eiXi > 0, e^x^ > (2.34) 

where 6j = ±, % = 1, . . . , N. Then, for x±, . . . , x^ > and xi, . . . , xn 2 by 2 different, we define 

$(xi,...,xjv) = <P ei '" eiT (eiXi,...,e N x N ) e ei ® • • • <g> e tN (2.35) 

ei,...,ejv=± 

where e + = ( _ J and e_ = ( ? ) . Here we stress that the original wavefunction y> is defined 



v o j — y i / 

for both signs of the Xj's and that the wavefunction $ contains the same physical information 
but is defined for x%, ... ,x^ > only. The advantage of the latter is that it allows to impose 
all the conditions on the wavefunction (interactions between particles in the bulk and boundary 
conditions at the impurity) in a very compact form as we shall see. 

Given a tensor product of spaces, (C 2 )®^, we define the action of a matrix M G End(C 2 ) on 
the k-th space by 

M [k] =1® •••®I®M®I® (2.36) 

fc-l N-k 

Therefore, the impurity conditions are, 

$( Xl ,...,x N ) = -i{U + I)W d th ${x 1 ,...,x N ) for x k ^0+ , l^k^N(2.37) 
The natural generalization of the bulk conditions read, for Q G &n and 1 ^ % ^ JV — 1, 

• • , Ql =4 (l+1) =5^* ^ $ (^ • • ■ XQi = XQ{i+1) (2-38) 

= Qf t Q- 1 [(d XQi - d XQ(i+1) ) + 2g] $(x 1; . . . ,a*)| XQl=x - (i+i) (2-39) 



i j ■ 



Q is the usual representation of the element Q G 6jy on (C )® . Namely, denoting by E. 
i,j = 1,2 the matrices with 1 at position and elsewhere, one has 

(2.40) 




kl=l 



Then using T{Tj = T(Tj and (|1.9p . it is easy to get Q for any Q G &n since an arbitrary permutation 
can always be decomposed in transpositions. At this stage, we have explicitly formulated the N- 
body problem corresponding to the combination of the two systems as described in Section [TJ 
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Let us make the ansatz for $: in the region < Xq\ < • • • < xqn, Q £ &n, it is represented 



by 



$q(xi, . . . , Xat) = ^ exp[^(/cpiXQi H h fcpTVXQTv)] QAp(Q) 

Pe 2Bjv 



where 



ei,...,ejv=± 



(2.41) 



(2.42) 



Again, the eigenvalue problem is simply solved by E = YliLi k h 
Inserting in (12. 3 7p . one gets 

A PRl (Q) = Z [1] (-kpi) A P {Q) 
where Z is given by (12. 61) . From relations (I2.38f) . (12.391) . we get for 1 ^ j ^ N — 1 



AptAQ) 



k Pj - kp(j +1 ) + ig 



((k Pj - kp {j+1) )Ap(QTj) - igA P (Q)) 



(2.43) 



(2.44) 



To get an analog of (I2.29p . we introduce an ordering on & N by associating to each element Q e & N 
an integer [Q] e {1, . . . , AH} so that Q be the [Q] th element of the ordering list. Next, we define 



(2.45) 



where e 



IQ] 









\[Q] 







1 













w 





e C m so that A P (Q) is just (A P ) 



[QY 



Thus, the relations (12.431) and (I2.44p take the compact form 

Apr^ = Zi(—kpi) Ap 

and for 1 ^ j < N - 1 



A 



PT, 



Yj(k 



P j — kp(j + i))Ap 



(2.46) 



(2.47) 



where the matrix elements of Z and Yj read (recall that these matrix elements are matrices 
themselves acting on (C 2 )® N ) 



Z i( k )[Q),lQ'} = zW ( k ) S [QUQ'] 

Y j( k )[QUQ>] = ~j—- g ( k 5 \QTMQ'\ ~ z 9 6 [QUQ']) f ® N 



(2.48) 
(2.49) 
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Since our construction is based on 2Utv, the Bethe ansatz solution is consistent if Z\ and Yj satisfy 



Yj(k)Yj(—k) = I m ® I® N , Z x {k)Z x {-k) = l m ® I 0iV (2.50) 

Yji^Yj+xih + k 2 )Y 3 {k 2 ) = Yj+xi^Yjih + k 2 )Y j+1 (h) (2.51) 

Y 1 (k 1 - k 2 )Z 1 (k 1 )Y 1 (k 1 + k 2 )Z x {k 2 ) = Z 1 (k 2 )Y 1 (k 1 + k^kjYrikt - k 2 ) (2.52) 

Y j (k 1 )Y e (k 2 ) = Y e (k 2 )Y j (k 1 ) for \j - £\ > 1 , (2.53) 

Z^OY^) = F J (A; 2 )Z 1 (A; 1 ) for j > 1 (2.54) 



where Ijyj the iV! x AT! unit matrix. Relations (I2.50p are usually called unitarity conditions while 
(I2.5ip is the celebrated quantum Yang-Baxter equation [7J [8]. Relation (12.521) is again our gen- 
eralized reflection equation. One can check that these relations hold true by direct computation 
(whatever the values of g and £, [i and v defined in (11. 6p ). finishing our argument about the exact 
solvability of our Y-particle system. Starting from Aid and using (12.461) and (12.471) repeatedly 
one gets the eigenfunction. 

3 Bethe ansatz: spectrum in the presence of an impurity 

In the previous section, we showed that the energy problem reads 

N 



where the fc's are the momenta of the particles. It is known that the complete use of the Bethe 
ansatz entails that the fc's are the solutions of the so-called Bethe ansatz equations. From these 
equations, it is possible to get some insight in the energy spectrum of the problem. The usual 
approach is to enclose the system in a finite region of space. One can imagine two types of 
conditions at the border of the finite region. In one dimension, one can put the N particles on a 
circle requiring periodic (or even anti-periodic) condition. This was the choice made in [IJ where 
the properties on the whole line were subsequently extracted through the so-called thermodynamic 
limit. An alternative approach is to enclose the particles in a box requiring the vanishing of the 
wave function on the walls of the box. This was explore e.g. in [9]. 

3.1 Bethe ansatz equations for particles on a circle 

Let us imagine that the N particles live on the interval [-L, L] centered for convenience on the 
impurity. In terms of the original wavefunction ip, the periodic (resp. anti-periodic) condition on 
the £-th particle, 1 ^ £ ^ N, reads, 





ip(x!, . . .,x ( -i,L,x e+1 . . .,x N ) 

(f (Xi, . . . , 2^-1, L, X£ + i . . . , Xjv) = 



= 9 ip(xi, Xi- t , -L, x e+1 x N ) 
8 (p'ixt, xt-x, -L, x i+1 x N ) , 




(3.3) 
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with 6 = 1 (resp. 6 = — 1). Introducing 

* = «(!J). < 3 - 4 > 

and using the tensor notations (I2.36p . equations (I3.2p and (13.31) can be equivalently written in 
terms of $ as 

. . . , xt-i, L, x e+1 ...,x N )= a [e] . . . , x e -i, L, x e+1 ...,x N ) (3.5) 

d xt $(xi, x e -!, L, x e+ i . . . , x N ) = -a [e] d xi $(xi, x £ -i, L, x e+1 . . . , x N ) (3.6) 

Invoking the Bethe ansatz solution (12.411) for some Qt e &n such that Qg(N) = £, one gets (noting 
that Q t <rM = a [N] Qfi 

e lk ™ L A P (Q e ) + e- lk ™ L A P R N (Qi) = <r [N] (e ikpNL A P (Q £ ) + e~ ik ™ L A PRN {Q t )) (3.7) 

e* k ™ L A P (Q e ) - e~ ikpNL Apr n (Qg) = -a™ {e ikpNL A P (Q £ ) - e~ ikpNL A PRn (Qi)) (3.8) 

This entails 

e 2lkpNL A P (Q £ )-a^A PRN (Qi)=0 , e=l,...,N (3.9) 
that is, in terms of A P as defined in (12.451) 



e 



2ik PN L Ap _ zA PRn = , E = Ijvi <g> a [N] . (3.10) 



This holds for any P G 2Uat yielding a priori 2 N N\ different equations. In fact, let us show that 
we only need to consider N of them by proving that if (I3.10p holds for A P then it holds for Apt > 
j = 1, . . . , N - 2, A PRl and A PRn . For j = 1, . . . , N - 2 

A PTj = Yj(k Pj - k P r j+1 ))A P (3.11) 



e 



-2ikpiyL 



<Yj{k Pj - k P(j+1) )XAp Rlt (3.12) 

= e-^^S^fe-fcp^))^ (3.13) 

= e-^^EAp^. (3.14) 

= e~ 2ikpNL J:A PTjRN (3.15) 

where we used Yj(k)Tl = EYj(fc) and -RatI} = TjRn. The proof for the other two cases is similar 
and requires Z 1 (k)T l = ZZi(k), i?7v#i = #i#tv and E 2 = I M ® I 0Ar . Since 7}, j = 1, . . . , N - 2 
and i?i are the generators of Wjv-i of cardinal 2 N ~ l (N — 1)!, adding i?^ brings the number of 
elements to 2 N (N — 1)!. Therefore, quotienting Wn by this set, we are left with N different 
elements: Sn = Id and Sj = Tj . . . T N _i for j = 1, . . . , N — 1. 
Now using (12.46P and (12.471) repeatedly, one has 

= yjv_i(-A:jv - • • • ii(-A;i - 

x ^(-fy)*!^! - . . . Y N ^{k N - kj)A Sj , (3.16) 
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and 

A Sj = Ytf-^kj -k N )... Yjikj - k j+ i)A Id . (3.17) 
Let us introduce the matrices TZj for j — 1, . . . , N as 

TZj = Yj(k j+1 - kj) . . . Y N _x(k N - kj) >: Kjv_i(-A;jv - %) • • • - kj) 

x Y}_!(— - kj) . . . Y x {-kx - k j )Z 1 {-k j )Y 1 {k l - kj) . . . ^-1(^-1 - kj) (3.18) 

Applying all these results in (13.101) . we are now in position to state the main result of this paper: 

Proposition 3.1 The wavefunction of our exactly solvable model is completely determined for a 
given vector Aid, using relations (STJQ) and \2.J^J\) to find Ap for any P G 2Hat. In turn, Aid is 



the common eigenvector of the matrices TZj with the eigenvalues e j respectively, j = 1, . . . , N: 

TZj A Id = e 2lk > L A M . (3.19) 

This entails in particular the following constraints 

det [TZj - e 2ik > L l m <g> l m ] = , j = l,...,N. (3.20) 

These are the impurity Bethe ansatz equations constraining the allowed values of the momenta of 
the particles. The presence of Z in TZj accounts for the effect of the impurity on the dynamics of 
the system while the matrices Yj contain the interaction effects. 

The proof goes as follows. First, relation ( 13. 19ft is a direct consequence of ( 13.101) . ( 13. 16ft and ( 13. 17ft . 
The fact that Aid is the common eigenvector of all these matrices follows from 

TZjTZ e = TZ e TZj , j,£ = l,...,N. (3.21) 

This equality, albeit tedious to establish, holds thanks to relations (12.50p - (12.54j) together with 

[Yj{h), Yj{k 2 )} = , for all k u k 2 . (3.22) 

3.2 Fixing the statistics 

So far, we said nothing about the statistics of the particles under considerations (on purpose). 
Indeed, our setup can be accommodated along the lines of [7j to allow for arbitrary statistics. Here, 
to get more insight when dealing with the Bethe ansatz equations, let us choose the statistics of our 
model. For bosons (resp. fermions), the wavefunction should be symmetric (resp. antisymmetric) 
under the exchange of any two particles. In terms of $, this reads, for 1 ^ % < j ^ N, 

$(xi, . . . , x h . . . ,Xj, . . . , x N ) = r $(xi, . . . , Xj, . . . ,Xi, . . . , x N ) (3.23) 

with r = +1 for bosons and r = — 1 for fermions. In turn, this yields an additional relation 
between the coefficients Ap(QTi) and Ap(Q): 

A P (QTi) = tA p (Q) . (3.24) 
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As a consequence, all the matrices Yj(u) become proportional to the identity, the multiplication 
factor being y T (k) given by 



rk — ig 



(3.25) 



k + ig 

When r = — 1, we recover that the 6 interaction is trivial for spinless fermions. The N impurity 
Bethe equations are equivalent to 



det 



II V T ( k 3 + k m )y T (k, - k m ) exppikjL) - Z^(-k 3 )a^ 



, l^j^N (3.26) 



Let us denote Z\(k), z 2 (k) the eigenvalues of Z(k) then, for N > 2, Z^(k)a^ has four different 
eigenvalues z\{k), —Zi(k), z 2 (k) and — z 2 (k), each of which is 2 Af_2 -fold degenerate. For N > 2, 
the N equations (I3.26f) are in turn equivalent to 



exp(2%L) Yl V T (kj + k m )y T {kj - k m ) = Xj(kj) , 1 < j < JV 



(3.27) 



where for each j, Xj(kj) takes one of the four possible values Z\{— kj), —Zi(-kj), z 2 {—kj) or 
—z 2 (—kj). Then, to find the complete spectrum, we need to solve the set of N equations (13.271) 
for the iV unknowns {kj}. Since, for each equation, four choices of Xj are possible, there are 4^ 
different sets of iV equations to solve. 
For N — 1, the equations read 



exp(2ikL) = si(—k) or exp(2ikL) = s 2 {—k) 



(3.28) 



where Si(fc), s 2 (k) are the eigenvalues of aZ{k). In the previous equations, we isolated on the 
left-hand side the usual terms corresponding to the interaction and on the right-hand side the new 
part arising from the presence of the reflecting and transmitting impurity. 



3.3 Bethe ansatz equations for particles in a box 

In this paragraph, instead of putting the particle on a circle, we are going to let them live in a 
box [-L, L\. Thus, we impose the following conditions for 1 ^ i ^ N 

$(xi, . . . , x e -i, L, x £+1 x N ) = , (3.29) 

and we follow the above analysis along the same lines. The linear system of equations relating the 
2^^! coefficients in each As d , j = 1, • • • , N, take the form 

KjAid = -e 2ik > L A Id (3.30) 

with IZj as defined in (13. 18j) replacing S by the identity here. In this case, the impurity Bethe 
equations read 

det [Kj + exp(2ikjL) l m ® l m ] =0 , j = l,...,N. (3.31) 
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These are the Bethe ansatz equations for our system with the particular box conditions under 
consideration. If we fix the statistics for bosons or fermions as in the previous paragraph, these 
equations take the simpler form 

exp(2z%L) ]JyT(k j + k m )y T (k j -k m ) = j j (k j ) , 1 ^ j < N (3.32) 

where for each j, jj(kj) takes two possible values : z\(— kj) or z 2 {— kj), yielding 2 N different sets 
of N equations. This setup will be used in the examples to illustrate the case of a parity-breaking 
impurity. 



4 Detailed study of selected examples 
4.1 Recovering previous results 

Let us show how to recover the historical results of [H] which are the analog of the results of [T] 
when the particles are confined in a box. 

We recall briefly the setup of [9] and adapt ours to reproduce it. M. Gaudin considered a 
gas of N bosons with 5 interaction enclosed in a box of length L. The symmetric wavefunction 
0(xi, . . . , xn) is required to satisfy 

4>(x! = 0,x 2 , ■ ■ ■ ,x N ) = 0, (4.1) 
4>(xi,x 2 , . • • ,x N = L) = 0, (4.2) 

in the region < x\ < x 2 < . . . < xn < L. 

Thus we must take (13.251) with r = 1. Then, the most natural idea that comes to mind 
to recover this setup from ours is to "fold" our system which lives on [— L, L] and tune the 
parameters of the impurity so as to make it a wall of the box at the origin. This goes as follows: 
for j — 1, . . . , N, we require 

<p(xx,...,x j ,...,x N )=<p(xi,...,-x j ,...,x N ) , < Xj < L (4.3) 

This global property has a direct consequence on $ as defined in (12.351) 

/ 1 \ 



. . . , x N ) = <f + - + (xi, . . . , X N ) 



)iV 



(4.4) 



V 1 / 



In other words, the representation is completely reducible and the only relevant wavefunction is 
ip + - + (xi, . . . , xn) which we identify to 4>(xx, . . . , xn)- The reducibility of the problem will show 
up again consistently in the rest of this paragraph (e.g. for Z(k)). It has to be related to the fact 
that we break the chirality of the impurity when we require (14.3 j) . that is we restore the parity 
invariance and the need for a two-dimensional representation disappears. 
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In this respect, we also expect (14. 3 p not to be compatible with the impurity conditions (12.371) 
in general. In fact, it is easy to see that the coefficients //, v in (I1.6P must satisfy 



^ _ n* = o , v + v* = Q. (4.5) 

When translated in terms of a,a,b,c,d in (12. 9p . this is perfectly consistent with the well-known 
characterization of a parity- invariant impurity i.e. a = d and a 2 = 1. 

Now to reproduce (14. ip . one just has to choose /j, — — 1 and v — 0. In particular, this gives 
= —I showing again the reducibility of the problem to a scalar representation where the 
impurity is to be seen as a purely reflecting wall with reflection coefficient equal to —1. 

Finally, taking 9 = — 1 in ( 13. 5\) . ( 13.6P yields ( 14.21) with no condition for the derivative as 
required. Again, the representation using a is reducible and one can see that the relevant eigenvalue 
for sigma is —1. Collecting all these settings, we end up with the following Bethe equations 

e 2ikjL = -Q A; m - kj - ig k m + kj - ig ^ j = i ^ N ( 4 6 ) 

, . k m kj -\- in k m -4- kj ~\~ io 

which are precisely those obtained by Gaudin in [9]. 

4.2 One and two particles with 5 impurity 

As a first step towards the understanding of the properties of our system, we pay special attention 
to the special cases of one and two particles in the presence of the well-known 5 impurity. The one- 
particle case is presented as a reminder and already displays the interesting features of degeneracies 
and bound states. Then, the two-particle exhibits the new properties arising from the impurity 
Bethe equations. 

The 5 impurity is characterized by 

tt r-^\ ( cos(0 isin(0 \ . r . v t . 

E7j=-exp(z0 • .% » ^e[0,7r) (4.7) 



zsin(^) cos(^) 

and the impurity Bethe equations constraining the momentum of the particle read 

exp(2i£;L) = 9 or exp(2ikL) = 9 kt ^^ + % (4.8) 

k tan £ — i 

The first equation reproduces the usual integer quantization of the momentum while the second 
shows how this quantization is controlled by £. 

Figure [1] shows the energy spectruuQ as a function of the tunable impurity parameter £. 

The constant energy levels correspond to the first equation and do not depend on the impurity 
parameter £ as expected. The other energy levels show the effect of the impurity on the spectrum. 
Of special interest is the lowest energy level which exhibits a bound state for £ > tt/2. This is 
consistent with the fact that the coupling constant to the impurity r), given by rj = 1/ tan£, be- 
comes negative. We have also plotted in Figure [2] the corresponding densities for different regimes. 



J all the figures are plotted in units of H 2 /2m, for g — 1 and for a unit length, L = 1, using Maple. 
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-1 1 

x 

Figure 1: Lowest energy level in terms of £ Figure 2: Density (=|<^(a;)| 2 ) for various val- 
for S impurity and 8 = 1 ues of £. 

The thick curve corresponds to £ = 0: the impurity is completely reflecting and no transmission 
occurs. For £ = n/3, the double-solid curve shows reflection and transmission for a repulsive 
impurity. The constant curve for £ = 7r/2 is very special since for this value of £, the impurity 
becomes trivial in the sense that the reflection vanishes and the transmission is just 1. This corre- 
sponds to the zero energy state. Finally, the thin curve represents the profile for £ = ll7r/12: this 
is the bound state whose profile gets sharper and sharper as £ — > n (infinitely attractive impurity). 

Let us move on to the case of two particles. The impurity Bethe equations read 

where each of the eigenvalues Xi(k), A 2 (A;) can be either ±1 or ^ ^ng-r ^ e display on Figure [3] 
the corresponding lowest energy levels. When the two eigenvalues are ±1, we obtain the constant 
energy levels. Otherwise, when one at least of the eigenvalues is ± fc ^|^ , the energy levels are 
decreasing functions of £. Again, for special values of £ (£ = 0, £ = vr/2), there are degeneracies 
which are lifted when we tune the impurity. Finally, for £ > it/2, the lowest energy levels give rise 
to bound states with the impurity (we recall that £ — ► tt corresponds to r] — > — oo i.e. an infinitely 
negative coupling constant). 

4.3 Asymmetric impurity in a box 

In this paragraph, we give an example of an asymmetric impurity i.e. an impurity which breaks 
parity invariance. We simply present the one-particle case for the box boundary conditions of 
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Figure 3: Lowest energy level in terms of £ for 2 particles and 5 impurity 



section 13.31 For convenience, we use the parametrization (12. 7p and the impurity is characterized 
by a single parameter as follows 

a — 1 , a = sin 2 w , b=— cosw , c = cosw , d = 1 , wG[0,7r). (4-10) 

Then the Bethe equations take the form 

2ikL (A; 2 — 1) cos w + iky/ 4 + cos 4 w 2ikL {k 2 — 1) cosw — ikyi + cos 4 w 

[k 2 + 1) cos w — ik(cos 2 w — 2) (k 2 + 1) cos w — ik(cos 2 w — 2) 

The two equations are never equivalent so that we do not observe level crossing (see Figure Hj). On 
Figure O the density profile for the lowest positive energy level shows the striking feature of this 
impurity for different values of w. One can clearly see the parity invariance breaking on the thin 
and double-solid curves (w = and Htt/12 respectively). Again for the particular value w = 7r/2 
(the thick curve), the impurity becomes reflectionless with trivial transmission equal to 1. Parity 
is restored and we observe a single excited mode in a box. 
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Figure 4: Lowest energy levels for a parity Figure 5: Density (=|<^(;c)| 2 ) for the first pos- 
breaking impurity itive energy level 

5 Conclusions and outlook 

In this paper, we presented and solved the one-dimensional problem of N interacting particles in 
the presence of an impurity. In the process of the Bethe ansatz for the wavefunction, doubling 
the dimension of the representation of the underlying Weyl group was a crucial ingredient with 
respect to previous approaches. This is reminiscent of the general RT algebras framework recently 
introduced and allows for an exact treatment of impurities. 

We also established the impurity Bethe equations controlling the energy spectrum. Although 
some basic understanding emerged from the study of the simple one and two-particle cases, a 
systematic study of the finite size Bethe equations as well as of the integral equations arising in 
the thermodynamic limit is needed. This will give important non-perturbative information on the 
effect of the impurity on the system. This should be compared to the perturbative and effective 
approach of Kane and Ficher [26] . A way to tackle this problem is to compute the sound velocity 
u and the Luttinger liquid parameter K directly from our impurity Bethe equations. The result 
will actually give the analogs of the renormalized u* and K* in Kane and Fisher's approach. This 
issue as well as other physical consequences deserve careful attention and will be investigated 
elsewhere. 

In any case, we already observed that a tunable impurity can lift degeneracies in the energy. It 
can also confine asymmetrically stationary states. In this respect, we emphasize that the present 
approach allows for the description of unusual asymmetric impurities (and not only the standard 
"delta impurity") whose effects for finite size systems and in the thermodynamic limit will also 
be addressed elsewhere. 
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